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Abstract 

We derive explicit formulas for probabilities of Brownian motion with jumps crossing linear 
or piecewise linear boundaries in any finite interval. We then use these formulas to approximate 
the boundary crossing probabilities for general nonlinear boundaries. The jump process can be 
any integer-valued process and jump sizes can have general distributions. Moreover, the jump 
sizes can be even correlated and/or non- identically distributed. The numerical algorithm is 
straightforward and easy to implement. Some numerical examples are presented. 
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1 Introduction 

Recently the first passage time (FPT) of jump-diffusion processes has drawn much attention in the 
literature, mainly due to its applications in finance and insurance. Some recent applications can be 
found in e.g. Cont and Tankov (2004), Kou and Wang (2004), Jiang and Pistorius (2008), Cai et al 
* Corresponding author 
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(2009) , Chiarella and Ziogas (2009), Kudryavtsev and Levendorskii (2009), Bakshi and Panayotov 

(2010) , Jeannin and Pistorisus (2010), Chi and Lin (2011), and Dong et al (2011). Given the 
important role of the FPT in many applications, however, the computation of the FPT densities for 
jump-diffusion processes turns out to be very challenging and so far very few exact solutions exist. 
In a remarkable work, Kou and Wang (2003) obtained the explicit form of the Laplace transform 
of the FPT for a Brownian motion (BM) with double-exponentially distributed jumps crossing a 
constant boundary. Dong et al (2011) obtained similar results for hyper-exponential jumps. Perry 
et al (2004) established the integral equation for the FPT density for BM with general jump sizes. 
In these proposed methods, the FPT density or probability have to be calculated through numerical 
Laplace inversion. More recently, Giraudo (2009) proposed an approximate solution for the FPT 
density for two-constant jump sizes. So far, no explicit form of the FPT distribution exists even 
for a constant boundary, let alone for general time-dependent boundaries. 

The calculation of the FPT distribution for diffusion processes is a long-standing problem. It is 
well-known that explicit formulas exist only for some special processes and boundaries. For general 
problems, there are mainly two analytic approaches. In the first approach, one establishes a certain 
integral or differential equation for the FPT density, and then solves this equation numerically to 
obtain an approximate solution. For example, the tangent approximation and other image methods 
are used by Strassen (1967), Daniels (1969, 1996), Ferebee (1982) and Lerche (1986), while a series 
expansion method is used by Durbin (1971, 1992), Ferebee (1983), Ricciardi et al. (1984), Giorno 
et al. (1989), and Sacerdote and Tomassetti (1996). 

In the second approach, one first establishes an explicit formula for piecewise linear boundary 
crossing probabilities (BCP), and then uses this formula to approximate the BCP for general nonlin- 
ear boundaries. This method was first proposed by Wang and Potzelberger (1997) for one-boundary 
problem and later extended to two-boundary problems by Novikov, Frishling and Kordzakhia (1999) 
and Potzelberger and Wang (2001). It has also been used by Novikov, Frishling and Kordzakhia 
(2003), Borovkov and Novikov (2005), and Downes and Borovkov (2008). One advantage of this 
approach is that the numerical computation is straightforward and the approximation error can 
be assessed and controlled. Moreover, this approach is easy to extend to deal with jump-diffusion 
processes, as is shown in this paper. 

Specifically, let (CI, (J r t)t>o,P) be a filtered probability space and (B t )t>o an J^-Brownian mo- 
tion. Further, let (N t ) t >o be a Poisson point process with rate A and (??i)igN a family of random 
variables representing jump sizes. Suppose that ((B t ), (N t ), (r^)) are mutually independent. Then 
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we consider the boundary crossing probability for the process 

X t = B t + ^2 Vi , (1.1) 

i=i 

where the second part of Xt is a so-called compound Poisson process. In this paper we first use 
the BCP approach of Wang and Potzelberger (1997) and Potzelberger and Wang (2001) to derive 
explicit formulas for the BCP of Xt crossing one- and two-sided linear boundaries respectively. Then 
we extend these formulas to piecewise linear boundaries and use them to approximate the BCP of 
Xt for general nonlinear boundaries. One advantage of this method is that it can also deal with 
more general jump component in Xt- For instance, Nt can be any integer-valued stochastic process 
starting with iVo = 0. Moreover, the jumps (^i)ieN axe not necessarily i.i.d. as long as the joint 
distribution of (rji, 772, • • • , %) is known for each k G N. This provides much more flexible models for 
real applications where the subsequent jumps are allowed to be dependent or even follow different 
distributions. Finally, although in this paper we deal with processes of the form (jl.ip explicitly, 
our method can be easily extended to processes of the form X t = fi(t) + aB t + ^=1 because 
the drift ii{t) and diffusion parameter a can be absorbed into the boundaries. 

This paper is organized as follows. In section 2 we derive explicit formulas for the BCP for one- 
and two-sided linear boundaries respectively. Then in section 3 we derive the formulas for piecewise 
linear boundaries and use them to approximate the BCP for general nonlinear boundaries. Finally, 
some numerical examples are given in section 4 and conclusions are in section 5. 

2 Explicit formulas for linear boundaries 

In this section we first consider the linear boundary crossing probability 

P(max{X, -as-b} < 0), (2.1) 

0<s<t 

where a and b > are constants. To simplify notation, throughout the paper we denote EL4] = 
E[1a] = P(-A) f° r anv measurable set A. Further, let 

P N = a(N s , < s < t), P B = a(B s , 0<s<t). 

A jump time is defined to be the time t > such that Xt 7^ X t -, where X t - = \\m s ^, t - X s . 
Subsequently, we denote a jump time by r and the corresponding jump height by r\ = X(t)—X(t~). 
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Thus, for each k € N, the first k jump times are 



r 1 =inf{t>0; N t ^N t -}, 
r 2 = inf{t >n; N t ^N t -}, 



T k = m£{t > r fe _ i; N t ^N t -}. 

Here similarly iV t - := lim s _j. t - iV s . For each A; £ N, let F k (ti, . . . , t k ) be the joint distribu- 
tion function of the jump times (ri,...,^) and G k {h\, . . . , h k ) be the joint distribution func- 
tion of the corresponding jump heights (r/i, 772, • • • , For each t > 0, the support of is 
S£ = {(si, . . . , Sfe); < Si < . . . < Sfc < t}. Further, denote the distribution of B t as 



X 2 



1 

$t(x) = — r= / e 2t dy. 



00 



Then we have the following results. 

Theorem 2.1 (One-sided linear BCP) For any constants a and b > 0, it holds 

00 „ k 

P( max {X s - as- 6} < 0) = V P(iV t = fc) / dF fc (ti, ti + t 2 , • • • , V <i) 

fc=0 ^fc i=l 

• / dG k (h 1 ---h k ) / ••• / H(l-e ) 

JR' J —00 J — 00 

• (Va^ + -|±L) - e-^ + i$(av^T- -^=)) d<M*i) . . .dS^x*), 
v V^+l v^+l 7 

where b\ = b, bi + \ = bi + ati — Xj — hi,i = 1, . . . , k; fa = minjatj + bi, ati + bi — hi}, i = 1, . . . , k; 

tk+i = t — Yli=i U > 0? anc ^ $0*0 := $10*0- Further, if N t is a Poisson process with rate A and (rji) 

are i.i.d. with distribution G, then 



(2.2) 



P( max \X S - as - b\ < 0) = V X k e~ xt dh dt 2 ■ ■ ■ / dt k 

°< s<t ^ Jo Jo Jo 

■ J^dG{h l )---dG(h k ) J ■J Xl{ l - e *i ) 

■ (* + ^g=) - e" 2 ^ 1 $ {a^i - -|g=) ) d* tl (xi) • • • d** (**)• 



(2.3) 
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Proof. First, since (Bt), (Nt) and (r^) are independent, we have 

P( max{X s - as - b] < 0) 

y 0<s<t ' 

= P ( max^ {B s + ^2 Vi - as - b } < 0^ 

OO n k n 

= J2?(N t = k) / dF fc (ti,...,J>) / dG(^i)...dG?(fc fc ) 

fc=0 5 fe i=l Rfc 



(2.4) 



•P( max (IL — as — 6) < 0, ... , max — as — 6+> /i;) < 

In order to calculate the term in the integration, we recall the following well-known results on 
Brownian motion (Siegmund, 1986, P. 375, G. Deelstra,1994): 

, f i I \ —2b(at-\-b — x) 

P( max {B s — as — bj > 0\Bt = x) = e * , for x < at + b, 
P( max {B s - as - b\ > 0) = 1 - $(ay/i + 4=) + e~ 2afe $(ay / t - -^=). 

Then we have 

fe 



P( max {B s — as — b] <0, . . . , max { B s — as — b+y hi] <0 



/ P( max { B s - as - b\ < 0\B t = xi) 
v o< s < tl l ' 



•P( max {B s -as-b+hi} < 0, max |i? s -as-6+y^ft-i) <0 i? tl = x x ) d$ tl (xi) 
v ti< s <ti+t 2 l ' T*u<8<t^ r-f 

z— 1 

/"ft 2b(at 1 +b- J1 ) / 

= / (l-e 'i )P max |S S - as - b 2 ] <0, . . . , 

J-oo \0<s<t 2 > 

k 

max {B s - as - b 2 + hA <0) d$ tl (xi). 
Iterating this calculation process yields 

fe 

PI max -f-B,— as — b\ < 0, ... , max {B s — as — b+y hA <0) 

= / •••/ IK 1 " 6 h )■ ^+7^^ 

J —OO J— OO -_i V 



1=1 



- e- 2 ^+^( av ^7 _ ^±L)) d* tl (xi) • • • d* tfc (x fe ). 
Substituting (|2.5p into (|2.4p we obtain the desired result (2.2). 



Remark 2.2 Formulas (12. 2p and (|2.3p are infinite sums due to the number of jumps of the process 
(J s ,s > 0) in interval [0,t). In practical calculations they have to be truncated at finite terms. 
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From (|2.4p it is easy to see that the truncation error at item k = n is bounded by Y^k= n ^(Nt = k), 
which is easily controlled by the behaviour of the process Nt. For example, if Nt is a Poisson 
process, then 



> P(N t = k)=y Vr-e < ) 

^ v ' ^ k\ ~ (n + 1)! 

k=n k=n 

Therefore, in practice one can determine a number n such that ( n+1 y is small enough, and then 
sum up the terms up to n — 1. Moreover, unlike in other methods such as by Kou and Wang 



(2003), controlling the truncation error in this way does not depend on the distribution of jump 
sizes r]i. In other words, the distribution of «j does not impact the accuracy of (|2.2|) and (|2.3[) in 
real calculation. 



Our second main result is based on the well-known formulas of Anderson (1960) [H (4.24) 
(4.32)], that for any t > and x < Sit + 71, 

72, S 2 ;t, x) :=P(3s < t,B s > 5 x s + 71, B u > 82U + 72, Vw € [0,s]|B 4 = x) 

00 

= ^ | e -(2A)[r 2 7i( 7 i+<5 1 t-x)+(r-l) 2 7 2 (72+<5 2 t-a:)-r(r-l){7 1 (72+(52t-a:)+72(7i+'5it-^)}] 



r=l 



-(2/t)[r 2 {7i(7i+i5it— x)+7 2 (72+i52t-a;)}— r(r— 1)71 (7 2 +<5 2 i— a:)— r(r— l)7 2 (7i+5ii— a;) 



(2.6) 



and 



X(7i>^i> 72,^2;*) := P(3s < t,B s > <5is + 7i,S„ > + Vu G [0,s]) 
'<Ji* + 71' 



1 - $ 



^ j e -2[r 7l -(r-l)7 2 ]Mi-(r-l)<5 2 ] $ ^ 
1 

2[r 2 (7i<5i+72<52)-)"(r-l)7i 5 2 -r(r— 1)<S 2 ] <j> ^ 



«5it + 2(r - 1)72 - (2r - l) 7l 



r=l 

e 



^it + 2r 72 -(2r-l)7i' 

Vt 



_ e — 2[(r— 1)71— r 7 2][() — l)<5i— r<5 2 ] 



1 - $ 



*it - 2r 72 + (2r - 1)71 



-2[r 2 (7i<5i+72<52)-»*(r-l)725i— r(r+ 1)71 5 2 ] 



1 - $ 



^it+(2r + l)7i -2r 72 

Vt 



}. (2.7) 



Then we have the following results. 
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Theorem 2.3 (Two-sided linear BCP) For any a, c G R and b, d > such that at + 6 > 
— (ct + d) for t > 0, it holds 

P(-(cs + d) < X s < as + 6,0 <s<t) 

= £p(jV t = A;) / dF fc (ti,--- ,X)*<) / dG(/n,...,%) 
fc=o ^ i=i ^ Rfc 

/•ft rPk * , v 

= / d$ tl (ari)--- / d$ tfc (ar fe ) fl -i?(6i,a,-c^,c|tj,Si) -i?(dj,c,-6i,-a|tj,-a;i)J 
■'ai -'"fe j = i 

• (l - x(Vn,o,-4+i,-c|< - - x(4+i,c, — fefc+i, — a|t - ( 2 - 8 ) 

i=i i=i 

w/iere 

6i = 6, = ati + bi - Xi - hi,i = 1, . . . ,k; 
j3i = mm{ati + fy, ati + 6, - /ij, i = 1, . . . , k; 
d\ = d; d i+ i = cti + di + Xi + hi, i = l,...,k; 
ai = minjctj + di, cti + di + hi}, i = l,...,k. 

Proof. Again by the independence of (B t ) t >o, {N t )t>o and (r/j)j e N, we have 

P(-(cs + d) < X s < as + 6,0 <s<t) 

N s 

= P(-(cs+d) < B s + ^r/i < as+b,0 <s<t) 
i=i 



= ^P(N t = k) / dF fc (ti,--- ,5^**) / dG(hi,...,/i fc )p(-(cs+d)<B a <os+6, 
fc=o i=i ^ Rfc 

<s<ti; . . . ;-(cs+d) < Bg + ^hi < as+b,^ti <s<tj. (2.9) 



i=l i=i 
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To evaluate the term in the above integral, we deduce inductively 

k 



P 



(-(cs+cf) < B s < as+b, <s<t\;... ;-(cs+d) < B s + '^hi < as+b, <s<tj 

i=i t=i 

- / P(-(cs+d) < B s < as+b, <s<ti\B tl = xi) 

J Ql 



• P( — (cs+d) < B s + h\ < as+b,ti <s<t x + 1% . . . ;-(cs+d) < B s 

k k 

+ ^2 h i< as+b^U <s<t\B tl = xi\d$ tl (xi)\ 



i=i i=i 



/ d$ tl (xi)--- / d$t k (xk) l[~P {- (cs+di) < B s < as+bi,0 <s<ti\B u = Xi ) 

k 

•P(-(cs+4+i) <B S < as+b k+ i,0 <s<t-^2u) 

i=i 

/ d$ tl (xi)--- / d^t k (xk)\^[i-^(bi,a,-di,c\ti,x i )-'&(di,c,-bi,-a\ti,-Xi) 

k k 

1 - x[h+l,a, -4+1, ~c\t -y~]tj) ~ x(4+i,c, -b k +i, -a\t - y]tj)\ (2-10) 

8=1 1=1 

where (|2.6p and f|2.T[) were used in the last step. The result then follows by substituting (|2.10p into 



3 Piecewise linear and nonlinear boundaries 

In this section, we first consider process Xt crossing a piecewise linear boundary. Suppose b(s) is 
a linear function on each subinterval of the partition = sq < si < • • • < s n —i < s n = t on [0, t]. 
Since the probability that any jump time falls on some points Sj is zero, without loss of generality 
we assume that any k > 1 realized jump times {uf, % = 1, 2, k} and points {sf, i = 1, n — 1} 
are distinct and therefore they form a partition < ti < ■ ■ ■ < t m = t, where m = k + n. Then 
analog to the proofs of Theorem 2.1 in section 2 and Theorem 1 of Wang and Potzelberger (1997) 
we can show the following result. 



Theorem 3.1 (Piecewise linear BCP) For the piecewise linear function b(t) defined on the par- 



Shao and Wang: Boundary crossing probabilities for Brownian motion with jumps 
tition = sq < si < • • • < s ra _i < s n = t, it holds 



P( max {X s - b(s)} < 0) = V] P(N t = k) [ dF k ( Ul , u 2 , . . . , u k ) [ dG k (h u h 2 , h k ) 

I ... j JJ (l - e - 2 ( 6 i-l- x i-l)( 6 i- x i)/(*i-*i-l)) 

. ~ Xm ) exp( - bm ^ brn ~ ~ 2&m + Xm ) ^ 

^ \/t ~ t'rn t — t m \ft t m / 

d$ tl (xi)d$ t2 _ tl (x 2 - ^l) • ■■d®t m -t m -i(Xm - Xrn-l), 

where hi = b(U) - Y^=i 1 {t i >u j }hj and fii = mm{bi, hi - hi}. 

Now we consider the BCP of process X t crossing a general boundary b(t). We first write 

P( max {X s - b(s)} < 0) = E[E( max {X a - b(s)} < 0|J^)], 

0<s<t 0<s<t 

where 

Pj = a(N s , < s < t,T k ,r](T k ),k> l). 

For any n > 1 and a partition < s\ < s 2 < • • • < s n = t, let b n (s) be the piecewise linear function 
connecting the points b(si),i = 1,2, ...,n and suppose 

max |6„(s) — b(s)\ — > 0, as n — > oo. 

0<s<t 

Then by the continuity property of probability measure, 

lim P(max{X s - b n (s)} < 0) = E[ lim Wmax{X 8 - 6 n (s)} < O^)] 

n-^oo 0<s<t ' n— >oo 0<s<t 

= E[E( max {X s - b(s)} < 0\Pj)} 

0<s<t 

= P(max{X s -fe( S )}<0). 

0<s<t 

Furthermore, by the results of Potzelberger and Wang (2001) and Borovkov and Novikov (2005), 
with probability one, we have 

P( max {X s - b n (s)} < 0\J*j) - P( max {X s - b(s)} < 0| Pj) = 0(\), 

o<s<t o<s<t n 

which implies that 

P( max {X 8 - b n (s)} < 0) - P(max {X s - b(s)} < 0) = 0(-U (3.2) 

0<«<t o<s<t n 

Hence (3.1) can be used to approximate the BCP for X t and general boundary b(t). 
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4 Numerical computation and examples 

In the previous sections, we provided various formulas for the BCP of jump-diffusion processes 
crossing linear or piecewise linear boundaries in the form of multiple integrals. Since these integrals 
are expectations of the integrants with respect to the probability distributions of jump process N t , rji 
and the multivariate normal variates B t . , they can be evaluated straightforwardly using the Monte 
Carlo integration method. As has been shown in Wang and Potzelberger (1997) and Potzelberger 
and Wang (2001), this method is very stable and produces fairly accurate results. In the following 
we illustrate this method using some examples. 

Specifically, we consider the Poisson jump process with rates A = 3 and A = 0.01 to approximate 
the no-jump scenario, and three different jump hight distributions: 

(1) double exponential (DE) distribution with density 

g(y) = p ■ me- my i {ym + (1 - p) ■ me my l {y<0} (4.1) 

and parameter values p = 0.5, rji = 1/0.10 and 772 = 1/0.15; 

(2) exponential (Exp) distribution with mean A = 0.15; and 

(3) Bernoulli (Ber) distribution 

g(y) = p ■ 1^=0.15} + (1 - P) ■ l{ w =-o.i5} ( 4 - 2 ) 

with p = 0.5. 

We calculate one-sided BCP for one constant (b(t) = 1), two linear (b(t) = ±0.5t+1.5) and three 
nonlinear boundaries (b(t) = 1 + 1 2 , a/1 + t, exp(— t)). For nonlinear boundaries, the approximating 
piecewise linear functions are based on n = 32 equally spaced points in the interval [0, t\. The cases 
of two linear boundaries are also considered by Kou and Wang (2003). 

We use the following procedure to compute the BCP: 1) generate the number of jumps Nt = k 
from the Poisson(At); 2) generate k uniformly distributed jump times u±, U2, in the interval 

[0,t] and the jump heights hi, h,2, hk from the above distribution G&; 3) generate the random 
sample from multivariate normal distribution of B tl , B t2 , B tm ; 4) evaluate the integrant in (3.1) 
say; 5) repeat the above steps TV times and calculate the Monte Carlo sample mean and standard 
deviation of the estimated integrant. 

In all cases, t = 1 and the Monte Carlo replications are iV = 200000. For comparison, we also 
include the case with A = 0. All numerical computations are carried out using computer package 
Matlab or R. R functions are available upon request. 
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The results for linear boundaries are given in Table 1, while those for nonlinear boundaries 
are in Table 2, where the Monte Carlo simulation standard errors are given in parentheses. The 
numerical results for constant and linear boundaries are consistent with those in the literature (Kou 
and Wang 2003). 

Table 1: BCP for linear boundaries with simulation standard errors in parentheses. 







A = 


A = 0.01 


A = 3 


b(t) = 1 


DE 


0.682826 


0.682841 


0.685138 






(0.000811) 


(0.000812) 


(0.000915) 




Exp 


0.682366 


0.679808 


0.106362 






(0.000812) 


(0.000816) 


(0.000639) 




Ber 


0.682884 


0.682599 


0.667722 






(0.000812) 


(0.000812) 


(0.000931) 


b(t) = 0.5* + 1.5 


DE 


0.941229 


0.942003 


0.938798 






(0.000404) 


(0.000401) 


(0.000466) 




Exp 


0.941859 


0.938322 


0.214821 






(0.000401) 


(0.000420) 


(0.000885) 




Ber 


0.941156 


0.941872 


0.933017 






(0.000405) 


(0.000402) 


(0.000487) 


b(t) = -0.54 + 1.5 


DE 


0.740624 


0.739806 


0.742868 






(0.000827) 


(0.000829) 


(0.000887) 




Exp 


0.738117 


0.735223 


0.120888 






(0.000830) 


(0.000836) 


(0.000689) 




Ber 


0.739535 


0.738522 


0.727180 






(0.000829) 


(0.000830) 


(0.000901) 



5 Conclusions 

We derived explicit formulas for probabilities of Brownian motion with jumps crossing linear or 
piecewise linear boundaries in any finite interval. These formulas can be used to approximate 
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Table 2: BCP for nonlinear boundaries with simulation standard errors given in parentheses. 









A = 


A = 0.01 


A = 3 


b(t) 


= l + t 2 


DE 


0.850851 


0.852495 


0.845398 








(0.000766) 


(0.000762) 


(0.000779) 






Exp 


0.851666 


0.847987 


0.167030 








(0.000763) 


(0.000772) 


(0.000823) 






Ber 


0.851581 


0.852696 


0.836923 








(0.000764) 


(0.000762) 


(0.000796) 


b(t) 


= y/l+t 


DE 


0.803963 


0.803150 


0.801449 








(0.000859) 


(0.000860) 


(0.000865) 






Exp 


0.804617 


0.797564 


0.141257 








(0.000857) 


(0.000870) 


(0.000768) 






Ber 


0.803450 


0.803659 


0.788627 








(0.000860) 


(0.000859) 


(0.000885) 


b(t) 


= exp(-t) 


DE 


0.439138 


0.439502 


0.454230 








(0.001079) 


(0.001079) 


(0.001084) 






Exp 


0.439860 


0.434360 


0.056572 








(0.001079) 


(0.001078) 


(0.000506) 






Ber 


0.437854 


0.437174 


0.430308 








(0.001078) 


(0.001078) 


(0.001078) 



the boundary crossing probabilities for jump-diffusions and general nonlinear boundaries. Unlike 
the FPT density method based on Laplace transform, our method does not reply on particular 
distributions for the jump components. In fact, the jump process can be any integer-valued process 
and jump sizes can have general distributions. Moreover, the jump sizes can be even correlated 
and/or non-identically distributed. Therefore this method is much more flexible and has much 
wider applications. The numerical calculation is straightforward and easy to implement. 

Finally we note that using the transformation method of Wang and Potzelberger (2007), it is 
possible to obtain the BCP for more general jump-diffusion processes of the form X t = Y t i Vii 
where Yt is a diffusion that can be expressed as a functional transformation of a BM. 
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